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ABSTRACT 

We have used GMOS-IFU and high resolution HST-ACS observations to map, in unprecedented 
detail, the gas velocity field and structure within the 0.7 kpc circumnuclear ring of the SBb 
LINER/Seyfert 1 galaxy NGC 1097. We find clear evidence of radial streaming motions associated 
with spiral structures leading to the unresolved (< 3.5 parsecs) nucleus, which we interpret as part of 
the fueling chain by which gas is transported to the nuclear starburst and supermassive black hole. 
Subject headings: Galaxies: active, Galaxies: kinematics and dynamics, Galaxies: nuclei 



1. INTRODUCTION 

Mechanisms of mass transfer in active galaxies from 
galactic scales down to nuclear scales, have been the sub- 
ject of many theoretical and observational studies (e.g., 
Shlosman, Begelman, & Frank 1990; Emsellem, Goud- 
frooij, & Ferruit 2003; Knapen 2005). Simulations have 
shown that non-axisymmetric potentials efficiently pro- 
mote gas inflow towards the inner kpc-scale regions (e.g., 
Engelmaier & Shlosman 2004). Further, recent observa- 
tions have revealed that structures such as small-scale 
disks or nuclear bars and associated spiral arms are fre- 
quently observed in the inner regions of active galaxies, 
suggesting they could be the means to transport gas from 
the kpc scale down to within a few tens of parsecs of the 
active nucleus (e.g., Knapen et al. 2000; Emsellem et al. 
2001; Maciejewski et al. 2002; Crenshaw, Kraemer, & 
Gabel 2003; Fathi et al. 2005). 

The nearby (14.5Mpc) LINER/Sey 1 galaxy NGC 1097 
is an important "laboratory" in this context. It hosts a 
large-scale bar as well as a secondary bar within its 0.7 
kpc circumnuclear ring (e.g., Shaw et al. 1993). Since the 
discovery of its double peaked Ha emission-lines (Storchi- 
Bergmann, Baldwin, & Wilson 1993), there have been 
several studies to follow the "fate" of the gas accumulated 
in the centre (e.g., Storchi-Bergmann et al. 2003). 

Recently Prieto, Maciejewski, & Reunanen (2005; 
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henceforth PMR05) have presented high resolution 
ground-based near-infrared images obtained with the 
Very Large Telescope. They reported the discovery of 
spiral structures within the inner 300 pc of the galaxy, 
which can be traced to within «10 pc of the nucleus (the 
limit of their spatial resolution). PMR05 argued that 
these spirals trace cool dust, and could be channels by 
which cold gas and dust are flowing to the nuclear super- 
massive black hole (hereafter SMBH). 

Here we present two-dimensional maps of the gas kine- 
matics and structure within the inner 0.5x1.0 kpc of 
NGC 1097, which appear to trace the fueling path to the 
unresolved nucleus. Our study is based on high signal-to- 
noise ratio spectra obtained with the Gemini Multi Ob- 
ject Spectrograph (GMOS-S) Integral Field Unit (IFU), 
which reveal the inner gas kinematics in unprecedented 
detail. We use these data in combination with high res- 
olution images obtained with the Advanced Camera for 
Surveys of the Hubble Space Telescope (HST-ACS) to 
cross-correlate the gas kinematics with the morphology 
of the nuclear spiral structure. 

2. OBSERVATIONS 

Our two-dimensional spectroscopic observations of 
NGC 1097 were obtained with the GMOS-IFU on the 
Gemini South Telescope, on January the 8th 2005, us- 
ing the R400_G5325 grating and r_G0326 filter (Gemini 
project GS-2004B-Q-25, PI: Storchi-Bergmann). These 
spectra cover the wavelength range 5600-7000 A at a 
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Fig. 1. — Sample GMOS spectrum of the nucleus illustrating the 
multi-Gaussian fit used to determine the velocities of the narrow 
lines and remove the underlying double-peaked broad H-alpha line. 
Each of the labelled narrow lines is modelled by a single Gaussian 
profile. The broad H-alpha line is represented by the 3 broader 
Gaussians. The individual Gaussian components are plotted below 
the spectrum, together with the residuals (dots) resulting from a 
simultaneous fit to the spectrum of all the Gaussians plus a fixed 
continuum level. The three "unmarked" Gaussians are a fit simply 
to remove nuisance structure due to the double-peaked broad Ha 
emission-line. 

spectral resolution of R«3500 (85 km s"^ FWHM). We 
observed three consecutive fields of angular dimensions 
5" X 7" each, covering in total 15" along the major axis 
of the galaxy (position angle PA=129°) and 7" along the 
minor axis, sampled by hexagonal lenslets with projected 
diameters of 0'.'2. Three exposures of 600 seconds were 
obtained for each field, slightly shifted along PA=40° in 
order to correct for detector defects. The spectra were 
reduced using standard tasks in Iraf. Products of the 
reductions are three wavelength and flux-calibrated data 
cubes, resampled to rectangular arrays of 50 x 70 spectra, 
each corresponding to an angular coverage of ff.'l x O'.'l. 

The HST-ACS image for NGC 1097 was obtained 
(GO 9782; PI: Axon) using the HRC mode (0'.'025 pix~^ 
and a field of view of 26" x 29"), through the FR656N 
filter. The ACS image was reduced and flux calibrated 
using the standard Calacs pipeline. 

3. RESULTS 

3.1. Spectroscopy 

We used the two-dimensional spectroscopic observa- 
tions to map the kinematics of the ionized gas in the 
central 0.5 x 1 kpc. Inspection of the individual spectra 
reveals the presence of the [O i] AA6300,6364, HaA6563, 
[Nii]AA6548,6583, and [S ii]AA6717,6731 emission-lines 
throughout most of the observed field. The velocity field 
was derived from the central wavelengths of Gaussians 
fitted to the narrow Ha component and to the [N ii] 
emission-lines, as these lines have good signal-to-noise 
ratio even at the borders of the fields. Within the inner 
0'.'5, we introduced 3 additional Gaussian components to 
fit and remove the underlying double-peaked broad com- 
ponent of the Ha line (see for example, Fig.^). However, 
this does not significantly affect the velocity field derived 
from the narrow lines. Since the velocity fields of Ha and 
[Nii]A6583 are identical within the uncertainties (« 10 
km s~^), we chose the latter to to represent the gaseous 
velocity field of the galaxy because of its higher signal- 
to-noise ratio. 



In the top panels of Fig. 13 we present the two- 
dimensional maps of the total flux in the bandpass, the 
integrated [Nii]A6583 flux, and velocity dispersion {a). 
The nucleus is at position (0, 0), while a part of the cir- 
cumnuclear 0.7 kpc star- forming ring is visible at the top 
of the panels. The cr map (corrected for the instrumental 
resolution) has circumnuclear values around 40 km s~ 
and decrease outwards to values of «15 km s~^, as one 
approaches the circumnuclear ring at about 7" south-east 
from the nucleus. 

The velocity field is shown in the bottom left panel of 
Fig. [5] It is clearly dominated by rotation with distor- 
tions indicative of significant non-circular motions. Simi- 
lar distortions have previously been observed in the large- 
scale gas kinematics of nearby disk galaxies (e.g. Visser 
1980; Zurita et al. 2004; Emsellem et al. 2006), where 
they were found to trace the galaxy spiral arms. In order 
to try to isolate these non-circular motions we have fitted 
an exponential thin disk model (Freeman 1970) to the ve- 
locity field using a method developed by Fathi (2004) . In 
the bottom middle panel of Fig. 13 we illustrate the best 
model, which has a deprojected maximum amplitude of 
«200 km s~^, a disk scale length of 3"5, a kinematic PA 
of 137° and disk inclination of 35°. This inclination is 
similar to that derived by Storchi-Bergmann et al. (2003) 
for the nuclear accretion disk, but is smaller than the 
value of 46° found from the large-scale gas kinematics 
by Storchi-Bergmann et al. (1996). The kinematic PA 
agrees with the value of 135° from Storchi-Bergmann et 
al. (1996) and Emsellem et al. (2001). We have tested 
the robustness of the model against extinction variations 
across the field by repeating the fit several times, mask- 
ing out the regions most affected by the distortions. We 
found no significant variations in the model parameters. 

We constructed a residual map by removing the disk 
model from the observed velocity field. The result is 
shown in the bottom right panel of Fig. |21 where it can 
be observed that the residuals display amplitudes of «50 
km s~^. Moreover, the velocity residuals are distributed 
in a spiral pattern, centered on the nucleus. We identify 
three main arms, designated 1, 2, and 3, which are out- 
lined by white dots superposed on the residual map. The 
residuals indicate redshifts along arm 2, but blueshifts 
along arms 1 and 3. As arms 1 and 2 are oriented approx- 
imately along the minor axis, and the south-west is the 
near side of the galaxy (Storchi-Bergmann et al. 1996), 
we conclude that the residuals indicate radial inflow of 
gas along these circumnuclear arms. This follows from 
the geometry of a configuration where the radial motions 
occur along the spiral arms which are in the plane of the 
disk. In this case, one expects the maximum line-of-sight 
velocities along the minor axis. Furthermore, this im- 
plies that the spiraling direction is clockwise. This simple 
setup does not fully explain the observed blueshift along 
the inner parts of arm 3, where outflows or individual 
gas clouds could complicate the observed velocities. 

3.2. Imaging 

We used two different techniques of image analysis to 
investigate the structures present in the image. The 
first is generating a "structure map" , which is based on 
the Richardson-Lucy image restoration technique (e.g., 
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Fig. 2. — GMOS-IFU data results: Reconstructed image formed by collapsing the 6300-6850 A wavelength range (top left); [Nll]A6583 
Flux distribution (top middle); together with the kinematic maps derived from the [Nll]A6583 emission-line, as well as the best fitting 
exponential disk velocity field model and residuals. The spiral features are delineated by white dots, with the numbers indicating the 
different arms as in Fig|3] In the bottom panels, red color indicates redshift and blue color, blueshift. All panels share the same orientation. 



Pogge & Martini 2002) . The second method we have ap- 
phed for enhancing the morphological structures makes 
use of the Iraf Stsdas/Ellipse ellipse fitting routine to 
construct a model for the galaxy. The structures are re- 
vealed in the residual image obtained by subtracting the 
model from the original image (see Fig.O. 

The original ACS image shows several dark filaments, 
in the form of spirals, which are enhanced in the pro- 
cessed images. A comparison with the J-band image of 
PMR05 (their Fig. 5), shows good agreement between 
the spiral structures seen in the optical and near-infrared, 
and supports their interpretation that the dark filaments 
map dust spirals which converge on the unresolved nu- 
cleus. The dark filaments form a pattern similar to that 
seen in the velocity residuals. However, a detailed com- 



parison of the velocity residual map with the ACS image 
shows that the spiral structure identified in the former 
is rotated relative to that traced by the dark filaments 
in the latter. To illustrate this, we have overplottcd the 
loci of the velocity residual spiral arms on the ACS im- 
age (Fig.|2Il. The radial streaming motions in the ionized 
gas occur between the spiral dust lanes. The latter are 
indicated by white crosses in Fig. |31 

4. CONCLUSIONS 

We have mapped the gas velocity field in the inner 
7"xl5"of NGC 1097 with GMOS-IFU at an unprece- 
dented spatial resolution. These data have enabled us 
to disentangle circular from non-circular motions in the 
ionized gas. The dominant motion is disk rotation, upon 
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Fig. 3.— The inner 10" X 10" of the HST-ACS image of NGC 1097 (left), the structure map (middle), and residual between the HST 
image and the elliptical model (right). The angular sampling of 0"025 pix~^ corresponds to f»1.8 pc pix~^ at the g alaxy. The dots delineate 
the spiral arms found in our kinematic analysis, labeled as in Fig.|2] which follow clockwise rotation (see Sec. 13.11 . The white crosses mark 
the dark regions from Fig. 5 of PMR05. 



which are superimposed velocity residuals that delineate 
three spiral arms converging onto the unresolved nucleus. 
We interpret the residual velocities - which reach «20% 
of the rotational velocity amplitude - as evidence for ra- 
dial inflow at ss50 km s~^ from the region interior to the 
circumnuclear 0.7 kpc ring, towards the nucleus. 

A similar spiral structure delineated by dust lanes is 
seen in both our optical HST-ACS image and the near- 
infrared images of PRM05. The radial streaming motions 
in the ionized gas occur between the spiral dust lanes. 
On larger scales, similar displacements between gas and 
dust are common in galactic spiral arms (Tilanus & Allen 
1991). This offset is attributed to loss of angular momen- 
tum of the gas as it passes through the spiral shock (e.g., 
Combes 1996). Assuming that the dust lanes trace the 
spiral shock, a similar mechanism may operate in the in- 
ner kpc of NGC 1097. As gas initially on circular orbits 
passes through the shock, some fraction of its kinetic en- 
ergy is thermalized, it loses angular momentum and falls 
toward the center of the gravitational potential. 

Taken at face value, the inferred streaming velocities 
would bring the gas from scales of 10 pc (the spatial 
resolution of our kinematic maps) down to the nuclear 
SMBH in approximately 200,000 yrs. Such a gas flow 
could also be responsible for triggering the recently dis- 
covered starburst activity in the immediate vicinity of 
the nucleus (Storchi-Bergmann et al. 2005). 

NGC 1097 seems to have all the necessary features to 



allow gas transfer from galactic scales down to the nu- 
cleus (e.g., Engelmaier & Shlosman 2004; Maciejewski 
2004). It has a large-scale bar, which is most likely the 
agent to bring gas down to the well known circumnuclear 
ring. It has a secondary bar, which further funnels gas 
toward the nuclear region, while here we have mapped 
streaming velocities along spiral structures, down to at 
least 10 pc from the nucleus. Further, the high resolution 
HST-ACS image suggests that these structures continue 
down to w3.5 pc scale. It is the first time that radial 
gas inflow has been mapped to such small radii from the 
central supermassive black hole. 
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